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Utilizing ASCA archival data of about 300 objects of elliptical galaxies, groups, and clusters of galaxies, 
^v^j ' we performed systematic measurements of the X-ray properties of hot gas in their systems, and compiled 

them in this paper. The steepness of the luminosity-temperature (LT) relation, Lx oc (kT) a , in the range 
of kT =1.5-15 keV is a = 3.17 ± 0.15, consistent with previous measurements. In the relation, we find two 
breaks at around ICM temperatures of 1 keV and 4 keV: a = 2.34 ±0.29 above 4 keV, 3.74 ±0.32 in 1.5-5 
keV, and 4.03 ± 1.07 below 1.5 keV. Such two breaks are also evident in the temperature and size relation. 
The steepness in the LT relation at kT > 4 keV is consistent with the scale-relation derived from the CDM 
model, indicating that the gravitational effect is dominant in richer clusters, while poorer clusters suffer 
non-gravity effects. The steep LT relation below 1 keV is almost attributed to X-ray faint systems of 
elliptical galaxies and galaxy groups. We found that the ICM mass within the scaling radius -R1500 follows 
the relation of M gas oc 7 l2 - 33 ± 07 from X-ray faint galaxies to rich clusters. Therefore, we speculate that 
"*^h \ even such X-ray faint systems contain a large-scale hot gas, which is too faint to detect. 

Key words: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — X-rays: galaxies: 
1 1 clusters — X-rays: ISM 
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, 1. Introduction the mass, M gas , of hot gas is more important and use- 

ful information for considering the characteristic of hot 
Elliptical galaxies, groups, and clusters of galaxies are gas. Arnaud and Raymond (1992) and David, Jones, 
known to contain a large amount of gravitationally bound and Forman (1995) derived the relation of the ICM mass, 
hot gas, which is bright in the X-ray band. Their gravi- M gas , and temperature, kT, although their sample was 



tational mass ranges over 3-4 orders from 10 12_13 M Q to small and the error of the ICM mass was large. A detailed 

10 16 Mq, and the temperature of the hot gas represents the relation of kT-M gas has recently been obtained by Mohr, 

depth of gravitational potential. Therefore, the tempera- Mathiesen, and Evrard (1999), who found a clear corre- 

ture dependences of the physical values of the hot gas are lation of M gas oc (kT) 1 - 98 within the radius -R500, within 

interesting in terms of hierarchical structure formations of which the averaged mass density is 500-times as high as 

galaxies and clusters and the origin of hot gas. the cosmic critical density. The information is still lim- 

Thc relationship between the X-ray luminosity, Lx, and itcd to rich clusters, and a relation for lower temperature 

the temperature, kT, of hot gas has been well investi- systems has not been obtained yet, because of their X-ray 

gated [e.g., Matsumoto et al. (1997) and O'Sullivan ct al. faintness. The ICM mass fraction is also an important 

(2001) for elliptical galaxies; Edge, Stewart (1991), David quantity to constrain the parameters of structure forma- 

et al. (1993), Ponman et al. (1996), Markevitch (1998), tion theories. However, the correlation of the ICM mass 

Allen, Fabian (1998), and Xue, Wu (2000) for clusters of with the cluster temperature is not straightforwardly pre- 

galaxies], since X-ray luminosities and temperatures can sented, because there is no clear integration radius that 

be easily measured. Another merit of this relation is a is well justified. Recently, radii within which the average 

weak dependence on the integration radius, and thus it mass density is higher by a certain factor than the cosmic 

can be well defined. critical density are frequently applied. This has theoreti- 

On the other hand, the X-ray luminosity of hot gas cal meaning since structure-formation theories show that 

depends on several properties, such as the hot gas tern- clusters of galaxies can form against the cosmic expansion 

perature, mass, and spatial distribution. Alternatively, within such defined radii; for example, Riso (Navarro et 
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al. 1995). However, the radius i?det, within which X-ray 
emission was detected, seems to be much smaller than 
i?i8o f° r elliptical galaxies and galaxy groups. Following 
Navarro, Frenk, and White (1995), Risq is ~ 1 Mpc for 
these objects although the radius i?det is at most 0.3-0.5 
Mpc for galaxy groups (Mulchaey et al. 1996a) and 0.05- 
0.2 Mpc for elliptical galaxies (Matsushita 2001). Since 
the radius i?dct depends on the sensitivity of instruments, 
systematic studies with the same instrument are necessary 
to discuss it. 

Concerning the origin of hot gas and the formation his- 
tory of systems, groups and poor clusters of galaxies are 
especially attractive because they are intermediate sys- 
tems between individual galaxies and clusters of galaxies, 
as indicated by the following three phenomena. First, 
the mass of hot gas in galaxy groups scatters widely re- 
gardless of their similar stellar mass, M sta r (Mulchaey et 
al. 1996a). The hot-gas-to-stellar mass ratio, M gas /M star , 
ranges from 0.01 to 5, and the lower values are similar to 
those of elliptical galaxies and the higher values arc typ- 
ical for clusters. Second, non-gravity heating is strongly 
suggested for these low-temperature systems. Ponman, 
Cannon, and Navarro (1999) discovered the entropy floor; 
the ICM entropy at the center of poor clusters exceeds 
the extrapolation from the relation between the tempera- 
ture and the entropy of rich clusters, suggesting substan- 
tial non-gravitational effects on the ICM in poor clusters. 
The entropy in the ICM is important to study the ther- 
modynamic history. A preheating model, which consid- 
ers significant energy input before accreting of gas into 
the cluster potential, has been suggested to account for 
this (e.g., Tozzi, Norman 2001), and claimed by some ob- 
servations (Finoguenov et al. 2001), although it is also 
inconsistent with the observations (Ponman et al. 2003; 
Mushotzky et al. 2003; Pratt, Arnaud 2003). The contri- 
bution of non-gravity heating is also consistent with the 
picture suggested by studies of the metal abundances; sig- 
nificant fractions of metals ejected from member galaxies 
have escaped from groups and poor clusters of galaxies 
(Renzini et al. 1993; Fukazawa et al. 1996; Fukazawa 
1997), indicating that galactic winds in the early galaxy 
formation epoch give a vast amount of energy to the sur- 
rounding hot gas. Third, Matsushita (1997, 2001) found 
that there are two types of elliptical galaxies: one exhibits 
a compact X-ray emission where the mass of hot gas is sev- 
eral percent of stellar mass, and the other shows an ex- 
tended X-ray emission that is represented by the double-/? 
model, and the mass of the hot gas is comparable to that 
of the galaxy groups. The latter type of elliptical galaxies 
has analogy with cD galaxy in groups, whereas the galaxy 
concentration around them is ambiguous. It is quite inter- 
esting why these two types exist. The above three issues 
are thought to have some relation with each other, and 
we can approach the unified picture by investigating the 
temperature dependence of the hot-gas properties from 
elliptical galaxies to rich clusters simultaneously. 

Here, we performed systematic measurements of the 
hot-gas properties as a function of the temperature, using 
ASCA (Tanaka et al. 1994) archival data without any se- 



lection criteria. Since observations of most objects were 
not performed by surveys, but proposed by many persons, 
the sample was not complete and was not easy to derive 
the luminosity, temperature, or mass function. The merits 
to utilize the ASCA data are as follows. The ASCA CIS 
is the most sensitive to the diffuse faint X-ray emission at 
the cluster periphery among the previous launched mis- 
sions, thanks to its stable low background, long exposure, 
and wide field of view. ASCA can measure a hard X-ray 
surface brightness that is less contaminated by the cen- 
tral cool component (e.g., Fabian et al. 1994). Moreover, 
the ASCA capability of resolving the emission lines leads 
to accurate measurements of emission integral of hot gas. 
Measurements with the same instrument and analysis pro- 
cedure are free from systematic calibration uncertainties 
among different instruments and analyses. This work is 
an extension of the ASCA results on 40 nearby clusters 
(Fukazawa 1997; Matsumoto et al. 2000). Throughout 
this paper, we assume the Hubble constant to be 50/i5o 
km s _1 Mpc -1 and qo = 0, and the errors represent the 
90% confidence range. 

2. Data Sample 

We utilized all of the ASCA data of elliptical galaxies 
and clusters. Several clusters were observed more than 
once, and we chose the observation with the longest ex- 
posure. The total number of objects that we identified 
as elliptical galaxies and clusters are 313, and are listed 
in table . Some of objects cannot be utilized for deriv- 
ing various correlations, due to an unknown redshift (17), 
insignificant detection (13), or contamination of the envi- 
ronmental X-ray emission (such as cluster emission around 
non-cD elliptical galaxies) (10). In the following, we cat- 
alogued the physical properties of hot gas for all of the 
objects that were significantly detected. 

Among 17 objects without an available redshift, the 
Fe-K line was clearly detected from 3 clusters, and we 
obtained the redshift by fitting the spectra with the red- 
shift parameter free. For the residual objects, we esti- 
mated them from the luminosity-temperature relation of 
Edge, Stewart (1991), L x = 1.0 x 10 43 (/cT/keV) 2 79 erg s" 1 
(2-10 keV). Here, the utilized temperature and the flux 
were derived by the following analysis procedure. The 
redshifts estimated here are summarized in table 3, as 
indicated by f, and were used to derive various physical 
quantities for these objects. There are 13 objects whose 
X-ray emission was not significantly detected, and we cal- 
culated the upper limits of the 0.5-2 keV flux within 6 
arcmin from each cluster center, assuming a power-law 
spectrum with a photon index of 1.5. In table 1, we 
summarized them. Among 10 objects that are contami- 
nated by the surrounding extended emission, we excluded 
NGC 4472, NGC 4406, NGC 4374, NGC 1404, NGC 499, 
and NGC 6034 in our correlation study, since the ambient 
cluster emission was too strong to constrain their X-ray 
surface brightness profiles. NGC 2865 was contaminated 
by the hard source, whose position was around 9 h 23 m 44 s 
and -23°08'50" (J2000). This object is extended in the 
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Chandra archival image, and its temperature and flux are 
4-5 keV and 1.7 x 1CT 12 erg s^ 1 cm~ 2 (0.5-10 keV), re- 
spectively. Considering the luminosity temperature rela- 
tion of clusters, it is thought to be an uncatalogued distant 
cluster whose redshift is 0.2-0.3. NGC 4291, CL 2236-04, 
and RX J1031. 6-2607 are also contaminated by the hard 
X-ray point source. Therefore, we excluded these four ob- 
jects in our study. For NGC 1316 and NGC 4649, the 
surrounding cluster emission is weak, and thus we ana- 
lyzed them, setting the cut-off radius by eye so that the 
surrounding emission would not affect the result. 

As a result, we analyzed 292 objects, among which ^50 
elliptical galaxies and galaxy groups were included. The 
redshift ranged from to 0.8, mainly 0-0.4, as shown in 
figure 1. Among them, we utilized 273 objects to derive 
various correlations. In figure 2, we show a plot of them on 
the Galactic coordinate. Most objects locate at the high 
Galactic latitude, and our sample lacks objects behind the 
Galactic plane. However, we expect that this selection 
effect did not affect our results. 

3. Analysis 

3.1. Data Reduction 

In the ASCA observations of the objects analyzed 
here, the GIS (Gas Imaging Spectrometer; Ohashi ct al. 
1996; Makishima ct al. 1996) data were all acquired in 
the normal PH mode, but the SIS (Solid-state Imaging 
Spectrometer) data were taken in various modes, such as 
the 4/2/1CCD Faint/Bright; the 4CCD mode was usual at 
the early ASCA phase, while 2CCD or 1CCD was frequent 
at the later phase because of CCD radiation degradation. 
Furthermore, the X-ray emission of objects is often over a 
small field of view of the SIS. For these reasons, we mainly 
analyzed the GIS data, and the SIS was utilized only to 
constrain the spectral parameters of the low-temperature 
objects, such as elliptical galaxies and groups of galaxies. 
The GIS data selections were performed on the condition 
of a minimum cut-off rigidity of 8 GeV c" 1 and a min- 
imum elevation angle of 5° above the earth rim. Flare 
events have been known to occur a few times per day in 
the GIS (Ishisaki 1996), and thus we excluded them. For 
the SIS, we further imposed the condition of an elevation 
angle larger than 25° above the day earth rim, and used 
events whose grade was 0, 2, 3, or 4. We used the SIS 
data of only the FAINT mode after 1995, 1996, and 1997 
for the 4, 2, and the 1CCD modes, respectively, and per- 
formed the RDD correction (Dotani ct al. 1996). In the 
analysis, we added all of the available data from different 
sensors, chips, modes, and pointings, separately for the 
GIS and the SIS after an appropriate gain correction. 

At first, we constructed GIS images of each observa- 
tion in the 0.6-7 keV, smoothed them with a Gaussian of 
a = 0'.5, and searched for contaminating sources within 25' 
of the detector center as follows. We searched pixels whose 
count rate was maximal, and marked it. We stamped it as 
an X-ray source. Next, we did the same thing, but by ex- 
cluding the regions within 3' of the pixels that had already 
been marked. If the marked pixel, except for the previous 
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Fig. 1. Redshift distribution of our sample objects. 
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Fig. 2. Distribution of our sample objects on the 
Galactic coordinate. 

one, was not within 4' of any marked pixels, we considered 
it as an additional X-ray source. We iterated the above 
procedure until the count rate of the newly marked pixel 
became less than the given value, which was the smaller 
one of either 2.5 x 10~ 5 c s _1 pix~ 2 or 3 c pix~ 2 ( 1 pix = 
0'.25 for thr GIS). We thus typically excluded 2'. 5 from 
X-ray source candidates, or even larger region for more 
extended or brighter sources for the GIS data. The ex- 
clusion of the same region was applied for the SIS data. 
After exclusion, we analyzed the GIS clean region on the 
detector coordinate, as defined in Ikebe (1995), which is 
defined so as to avoid a region where the background level 
is high or the calibration radio-isotope contaminates sig- 
nificantly. 

3.2. Background Subtraction 

In order to constrain diffuse X-ray emission, especially 
at the cluster periphery, we should subtract the GIS back- 
ground as accurately as possible. The background consists 
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Table 1. Lists of clusters with upper limits of flux in the observer frame at the 90% confidence level. 



Name F x (0.5-2 kcV) 

(erg s _1 cm~ 2 ) 



NGC 5018 


9 


.8 


X 


io- 


14 


GHO 1322+3114 


1 


3 


X 


io- 


13 


J 1888. 16CL 


5. 


9 


X 


io- 


14 


CL 0317+1521 


4 


5 


X 


io- 


14 


MS 1512.4+3647 


1 


.0 


X 


io- 


12 


PRG 38 


6 


9 


X 


io- 


14 


SCGG 205 


6 


9 


X 


io- 


14 


RGH 101 


9 


1 


X 


io- 


14 


3C 184 


8. 


5 


X 


io- 


14 


RX J1756. 5+6512 


1 


6 


X 


io- 


13 


3C 324 


5 


A 


X 


io- 


14 


PDCS 01 


2 


.8 


X 


10" 


14 


MS 0147.8-3941 


5 


.0 


X 


io- 


14 



of the cosmic X-ray background (CXB) and the intrinsic 
detector background (IDB). In practice, for both image 
and spectral analyses, we produced the background data 
set in the following way, which is almost the same as that 
described in Fukazawa et al. (2001), but a little different 
because the ICM emission is sometimes extended beyond 
the GIS full field of view. 

We first summed data of the ASCA Large Sky Surveys 
(Ueda ct al. 1999), conducted in 1993 December (AO-1 
phase) and 1994 June (AO-2 phase) over blank sky fields 
near the north ecliptic pole. Then, after Ikebe (1995), 
Ishisaki (1996), and Ueda (1996), we excluded regions in 
the GIS image where the count rate exceeds that of the 
surrounding region by > 2.5c. These regions correspond 
to faint sources with a 2-10 keV flux of > 8 x 10~ 14 erg 
s" 1 cm- 2 . The total exposure time of the GIS background 
data amounts to 233 ksec, ensuring negligible statistical 
errors. However, the derived background data cannot be 
used immediately, since the IDB level of the GIS has been 
gradually increasing by 2-3% per year, and it exhibits 
day-by-day fluctuations by 6-8% in the standard deviation 
(Ishisaki 1996; Ishisaki et al. 1997). 

Taking into account these effects, we estimated the IDB 
level of each pointing data individually, by assuming that 
the IDB spectrum and its radial profile arc both con- 
stant. Specifically, we derived four GIS spectra, denoted 
as S(E), B(E), N(E), and V(E); the former three were 
obtained from the on-source data, the blank-sky data, and 
the night-earth data, respectively. The last one is a sim- 
ulated spectrum of the ICM emission that was estimated 
by a full XRT+GIS simulation based on the ICM tem- 
perature and spatial distribution, both of which are first 
unknown. For a background estimation, accurate infor- 
mation of the ICM is not necessary. We then obtained it 
by ignoring the IDB time variation, by utilizing the GIS 
background data set described above, which is correct in 
a 0-th order approximation. A detailed description about 
how to derive the ICM temperature and spatial distribu- 
tion is given in the next subsection. 

These spectra were accumulated over the outer regions 



of the GIS field in the 6-10 keV energy range, to ensure 
that the CXB and V(E) are relatively minor compared 
to the IDB in S(E) and B{E). Note that N(E) con- 
sists solely of the IDB. We next fitted S(E) with a linear 
combination of B(E) + fN(E) + V(E), where / is a free 
parameter, and fN(E) represents the IDB difference be- 
tween the two epochs when S(E) and B(E) were acquired. 
Then, the fraction / turned out to be around 0.00-0.30 for 
any observations, as shown in figure 3. These values agree 
with the long-term increase of the GIS IDB (Ishisaki et 
al. 1997). By analyzing various ASCA data, we also con- 
firmed that this method can reproduce, within 5%, the 
GIS background spectra and radial profiles acquired at 
any time over 1993-1999. The background data obtained 
in this way were utilized in the subsequent image and 
spectral analyses. 

For SIS spectral analyses, we subtracted the back- 
ground in a conventional way, utilizing the archival SIS 
background set. 

3.3. Analysis Procedure 

As described in the previous subsection, the background 
estimation requires both spectral and imaging informa- 
tions of objects. Therefore, we first performed analyses 
by ignoring the gradual IDB increase to set / = 0, and 
obtained the spectral and imaging parameters somewhat 
roughly, but they were sufficiently accurate to constrain 
the IDB level. After estimating the IDB parameter, /, we 
again reanalyzed to constrain the spectral and imaging 
parameters more accurately. 

In the spectral and imaging analyses, we first deter- 
mined two radial extents for each object, i?iE6 and -R3E6, 
within which the GIS count rate was higher than 1 x 10 -6 
and 3 x 10" 6 c s _1 pix~ 2 , respectively, in the 0.9-7 keV. 
The GIS2+3 azimuthally averaged radial count rate pro- 
file after the background subtraction was used to deter- 
mine them. 
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Fig. 3. Distribution of the background correction factor, /, 
in the GIS data (sec text in detail). 

3.4- Spectral Fittings 

We here performed spectral fittings with the IDB factor 
/ free, and obtained the ICM temperature. The choice of 
integration radius for the spectral analysis followed the 
rule given in table 2, which was determined so as to avoid 
the central cool region. For objects whose spectrum is 
statistically poor, we set the inner radius to be smaller 
so as to include the bright center region. The ancillary 
response file (arf file) was calculated by averaging arf files 
at various positions, prepared by a step of 2'. 5 from the 
XRT optical axis. The spectral model applied here was 
the MEKAL model (Liedahl et al. 1995) multiplied by 
the photoelectric absorption. The column density was 
fixed to the Galactic value (Stark et al. 1992). In mea- 
surements of the temperature, the metal abundance was 
assumed to be 0.25 solar abundance, which refers to the 
solar photospheric values in Anders and Grevesse (1989), 
in order to minimize the error, especially for faint objects. 
Basically, we utilized only the GIS data, since the SIS cov- 
ers only the inner region of clusters due to its small field 
of view, or the SIS data at a later period often exhibit a 
severe RDD that changes a real spectral shape and cannot 
be completely corrected. On the other hand, for objects 
with kT < 2 keV, we performed simultaneous GIS+SIS 
spectral fittings because the SIS can well constrain the 
low-temperature component. In addition, we included a 
bremsstrahlung model for objects with kT < 1.5 keV, to 
reproduce the emission from low-mass X-ray binaries in 
elliptical galaxies (Matsushita et al. 1994), or the unre- 
solved excess hard X-ray emission from groups of galaxies 
(Fukazawa et al. 2001). Spectral fittings were performed 
in the 0.85-9.0 keV region for the GIS and 0.5-8.0 keV 
for the SIS. The thus-obtained results of the spectral fit- 
ting after considering the IDB increase are summarized 
in table . Besides 12 distant clusters and 2 poor clusters, 



the temperatures are well constrained with a < 20% ac- 
curacy. In figure 4, we compare the temperature of our 
results with that in Sanderson et al. (2003). Our temper- 
atures are systematically lower than those of Sanderson 
et al. (2003). Since Sanderson et al. (2003) determined 
the emission-weighted temperature by extrapolating the 
negative radial gradient of the temperature toward the 
central cooling region, the temperature is thought to be- 
come higher due to large emissivity at the center. 

Next, we performed a spectral fitting with the metal 
abundance free to constrain it. As a result, it was con- 
strained with a < 30% accuracy for about 27 objects. In 
figure 5 (left) , the obtained metal abundances are plotted 
against the temperature. Here, we do not plot objects 
whose error of the metal abundance is > 50%. Above 
kT > 2 keV, the metal abundance seems to negatively 
correlate with the temperature. This trend can be ex- 
plained by the effect of the high-mctallicity region at the 
bright cluster core, as claimed by Fukazawa et al. (2000), 
who showed that the Fe abundance at the outer region 
is almost constant to be 0.2-0.3 solar for any temper- 
ature (Fukazawa et al. 1998). Several reports on the 
metallicity increment at the center of clusters support 
this explanation (Ezawa et al. 1997; Fukazawa et al. 
2000; De Grandi et al 2001). In the choice of inner ra- 
dius in this analysis, we make it as small as possible 
so as to accumulate many photons. In order to look at 
this effect more effectively, we plot the emission-weighted 
metal abundance in figure 5 (right), which was obtained 
from the spectra over the cluster. The negative correla- 
tion of the metal abundance with the temperature be- 
comes prominent. The high-redshift clusters exhibit a 
similar trend, in good agreement with previous studies 
(Mushotzky, Loewcnstcin 1997; Matsumoto et al. 2000). 
Accordingly, when we study the metal abundance of dis- 
tant clusters that cannot be resolved spatially, we must 
pay attention to the effect of high metallicity at the cen- 
ter region. Below kT < 1 keV, the correlation disappears, 
and many objects exhibit a low metal abundance. This 
might be due to problems concerning the plasma model- 
ing and multi-temperature effects (Fukazawa et al. 1996; 
Buote 2000; Matsushita et al. 2003). 

3.5. Imaging Analysis 

After determining the ICM temperature, the az- 
imuthally averaged radial count rate profile was fitted with 
the singlc-/3 model. The column density and metal abun- 
dance were fixed to the same value as in the spectral anal- 
ysis of the temperature. We assumed isothcrmality, where 
the ICM temperature is that obtained by spectral fitting, 
and simulate the incident energy spectra with the plasma 
emission model by multiplying the photoelectric absorp- 
tion. Here, we applied the Raymond-Smith plasma model 
(Raymond, Smith 1977) rather than the MEKAL model in 
order to save time. The integration radii and energy band 
for imaging analyses followed the rule in table 2. The 
azimuthally averaged radial profiles of the X-ray counts 
were constructed by summing the GIS2 and GIS3 data. 
The center of the X-ray emission was set to the position 
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Table 2. Definition of radial and energy boundaries for spectral and imaging analyses. 



Radial boundary for spectral analysis 


Rsm > 10' 




3'-10' 


6'.75 < i?3E6 < = 


10' 




5' < R 3E6 <= 6' 


75 


9/ p/ 

z JX 3_E6 


R3E6 <= 5' 




l'-5' 


Boundary for imaging analysis 






Radius 


-R1E6 >— 15' 




0'-22'.5 


10' < i?iE6 < 15' 




0'-1.5r' 


-R1E6 <= 10' 




0'-10' 






Energy (kcV) 


kT <= 10(kcV) 




1.0 + 0.767 fcT 2 (lower energy limit) 






3.54 + 3.54 fcT( - k ° V - ) — - (higher energy limit) 


kT > 10(kcV) 




4.45-19.5 




2 5 
kT (keV) 

Fig. 4. Comparison of temperatures between our results and 
those of Sanderson et al. (2003). 

where the X-ray brightness became maximal. In order 
to model an X-ray surface brightness, we made use of a 
software system, called an ASCA simulator (Ikebe 1995), 
in which we produced the XRT+GIS angular response by 
using the Cyg X-l images actually observed with the GIS £ 
at 6 positions in the field of view (Takahashi ct al. 1995).° 
The actual radial profiles were compared with the sim- 
ulated ones through a \ 2 evaluation, where we assigned 
5% systematic errors to the background and the simulated 
profile, according to Ikebe (1995). 

We represented the ICM density distribution with a 
standard single-/? model (Cavaliere, Fusco-Femiano 1976) 
for all objects. However, an additional emission compo- 
nent is often required for poor clusters in the single (3 
model fitting. In such a case, we also applied a double-/? 
model that succeeded the representation of the X-ray sur- 
face brightness for some groups, rather than the simplc- 
[3 model (Ikebe et al. 1996; Mulchaey, Zabludoff 1998). 
The double-/? model consists of a small-scale central com- 
ponent and a large-scale component. The X-ray surface 





2 

kT (keV) 

Fig. 5. Fe abundance against the object temperature. The 
left panel plots the results for the spectral fitting in the default 
analysis region. The right panel plots for the whole cluster 
region. 
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brightness was expressed within a specified radius, i?ii m it, 
outside of which we assumed the ICM density to be zero. 
We determined i?n m it by increasing it until the \ 2 did not 
improve significantly. The surface brightness profile was 
thus calculated as 



^2 J ni{rfk{T)dl 



»=i 

. where r is the 3-dimcnsional radius, I is toward the line- 
of-sight direction, and A (T) is a cooling function calcu- 
lated with the Raymond-Smith model by assuming 0.25 
solar abundance; k is 1 or 2 for a single or double- (3 model, 
respectively, and i takes a value of 1 or 2 for a large-scale or 
small-scale component in the double-/? model. The square 
average of the electron density distribution, n|(r), is ex- 
pressed as 

91 -3ft 



1+ - 



Qi=l,2), 



is a central value and a,- is the core radius. For 



where ' 

th e doublc-/3 m odel, the electron density is expressed as 
\/n\(r) +n\(r). Since the spatial resolution of the XRT 
+ GIS is poor to constrain and fa of the small-scale 
component in the double-/? model, we fixed them to be 
a2 = kP c and fa = 0.7, which are averaged values 

of the small-scale components in Mulchaey and Zabludoff 
(1998). Therefore, the free parameters are i?ii m it, ^o,i, 
oi, and fa in the single-/? model fitting, and no, 2 becomes 
additional in the double-/? model fitting. 

Instead of building an automatic algorithm of % 2 mini- 
mization, we generated a model profile at each of the 8x8 
grid points in fa-a\ space, and calculated x 2 by adjust- 
ing the model normalization so as to archive the best fit 
between the data and the simulation. After this, we per- 
formed the same procedure by narrowing the fa—ai space 
around the above best-fit value. From the normalization, 
we calculated the central value, n 2 , j , by assuming the ratio 
of electrons to ions to be 1.18 and the luminosity distance 
to the object to be 



D 



zqo + (go - 1) (V2zq + 1-1) 
H q 2 



where go = 0.5 and Hq = 50h^o km s _1 Mpc -1 . For ellipti- 
cal galaxies and some galaxy groups, the hard component 
is significant in the spectra. Although we excluded the 
hard X-ray band in the image analysis, the contamination 
to the count rate must be taken into account. For this 
purpose, we estimated the photon count fraction of the 
thermal component and considered it in converting the 
surface brightness into the flux or n 2 , 1 . 

Examples of the best fittings are shown in figure 6, 
where the best-fit model is plotted as a solid line, su- 
perposed on the actual count-rate profile. The best-fit 
parameters and the lower limits of i?n m it are tabulated 
in table . The single-/? model fitting was successful in 
most cases with a reduced \ 2 of 0.8-1.2. The parame- 
ter errors for bright objects correspond to the parame- 
ter step, and thus the true errors may be smaller. The 



radius i?i; m it often corresponds to the radius of the GIS 
field of view of ~ 20'. fa and ai do not change beyond 
their parameter search step, as long as i?ii m ;t is within the 
confidence range. There are 57 objects whose X-ray sur- 
face brightness profile is well represented by the double-/? 
model, rather than the single-/3 model. In the following, 
we adopt the result of the double-/? model for these 57 
objects, which are summarized in table . 

From the tables, it is found that we can determine the 
core radius, oi, down to even ~ 1' with the ASCA data 
in this method. We compare the best-fit /? with that in 
Mohr, Mathicsen, and Evrard (1999) and Sanderson ct al. 
(2003) in figure 7. It can be seen that our results for /3 
arc in agreement with previous ones obtained in the soft 
X-ray band. However, our (3 and core radii are systemati- 
cally smaller than those of previous studies, by 5-10% and 
10-50%, respectively. Since these two quantities correlate 
with each other in the fitting, this could be an instrumen- 
tal artifact caused by, e.g., uncertainties in the XRT-PSF, 
differences in the fields of view, and so on. On the other 
hand, this might indicate that the X-ray surface bright- 
ness in the hard X-ray band is different from that in the 
soft X-ray band. A more detailed discussion on the X-ray 
surface brightness is beyond the scope of this paper. 



3.6. Temperature Determination under 
Consideration of the XRT-PSF Effect 



Full 



After obtaining imaging parameters, we further fitted 
the ICM temperature by considering the XRT scatter- 
ing effect, since the temperature obtained by simple spec- 
tral fittings might still be affected by XRT-PSF effects for 
high-tcmpcrature clusters (Takahashi ct al. 1995) that de- 
pend on the X-ray emission distribution. In this study, we 
excluded elliptical galaxies surrounded by the intraclutcr 
medium. Since the vignetting function strongly depends 
on the photon energy, we must deal with not the emission- 
weighted single spectrum, but the spectral sets consisting 
of spectra at several radii for determining the tempera- 
ture. We obtained two or three spectra at an appropriate 
annulus for both the data and simulation, and calculated 
the x 2 value of the difference between the data and the 
simulated spectral set by adjusting the normalization. In 
preparing the simulated data, we modeled the ICM den- 
sity distribution by the isothermal single-/? model with the 
obtained single-/? model parameters. We scanned the tem- 
perature from 0.4T S to 1.6T S by steps of 0.2T S , where T s 
is the temperature obtained in the simple spectral fitting, 
and calculated the chi-square value \ 2 f° r each tempera- 
ture. Then, the best-fit temperature, Tj, was determined 
by fitting the T-\ curve with the 2nd polynomial func- 
tion, as shown in figure 8. We compared the results with 
the temperature obtained by the simple spectral fitting, 
in figure 9. The temperature differs by at most 10% be- 
tween the two estimations; hereafter, we thus adopt the 
temperature obtained by simple spectral fittings. 

3. 7. Flux, Luminosity, and ICM Mass 

We derived an observed X-ray flux, Fx, and luminos- 
ity, Lx = ^irD 2 Fx, in the 0.5-10 keV energy range in the 
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Fig. 6. Examples of fittings of the radial count rate profile. 
The solid line represents the best-fit model and the cross 
points are the GfS data of the NGC i600 group. (Left) 
single-/? model, (right) double-/3 model. 
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Fig. 7. Comparison of f} of our results with those in Mohr, 
Mathiesen, and Evrard (1999) and Sanderson et al. (2003). 
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Fig. 8. Example of a \ 2 scan in the spectral fitting, 
considering the XRT-PSF effect for MKW 4. The ver- 
tical axis represents the \ 2 value. 
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Fig. 9. ComparisonSSsl|hcQcmfQrature obtained here 
T; m with that obtained in a simple spectral fitting (T sp ). 
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rest frame, by integrating the simulated emission spectra 
of the best-fit model within the best-fit i?ii m it- In order 
to consider the emission out of the GIS field of view for 
objects with i?ii m it > 18', we integrated the flux up to the 
i?^ it (fcT), which is an average detection radius for the 
specific temperature (figure 14): 0.2(/cT/lkeV) 1 - 70 Mpc 
for fcT<3keVand 1.3(fcT/3keV) 823 Mpc for fcT>3keV. 
This correction is typically at most 1.0-1.2, but ~2.0 for 
several nearby bright clusters, such as Centaurus, Perseus, 
and Virgo. We also calculated the bolometric luminos- 
ity by using the emissivity of the plasma model in the 
0.001-200 keV band in the rest frame. For the double-/? 
model, we derived it separately for the small and large- 
scale components. In order to estimate the error range of 
each quantity, we fitted the radial count rate profile by 
the (3 model for all of the imaging parameter sets within 
the 90% confidence range. We then calculated the flux 
and luminosity for each parameter set, and obtained the 
permitted range of each quantity. The obtained results 
of flux and luminosity, together with the central X-ray 
surface brightness, are given in table . The errors are typ- 
ically less than 15%, which are primarily attributed to the 
uncertainties of (3, the core radii (oi, 0,-2), and i?ii m it- In 
figure 10, we compare the bolometric luminosity obtained 
here with that in Xue and Wu (2000), who collected the 
available results of various authors. Both correlate well, 
but our results arc systematically lower by ^20%. One 
reason is that Xue and Wu (2000) corrected the observed 
luminosity by taking into account the undetected emis- 
sion at the periphery, while our results were calculated 
within the detection- limited radius. This effect explains 
the difference of 10-15%. 

We checked the consistency for the ICM mass between 
ours and the previous results of Mohr, Mathiesen, and 
Evrard (1999) and Sanderson et al. (2003). The integra- 
tion radius was chosen from this literature for the same 
objects analyzed here. Since the ICM mass, itself, is not 
available in Sanderson et al. (2003), we calculated it from 
the /3-model parameter in their paper. The ICM mass 
density was calculated from p sas (r) — fi e m p n e (r), where 
we assumed /i e = 1.167, based on the single-/? model fit- 
ting. In figure 11, we compare our results with those of 
two references. The systematic difference is seen against 
both references, although our results are in the middle of 
the difference between the two papers. 

The derivation of the ICM mass is not simple, because 
it strongly depends on the integration radius. The well- 
defined radius (Evrard 1997) r§oo, within which the aver- 
age mass density is 500-times as high as the critical den- 
sity of = 1-56 x 10~ 29 /i5Q g cm~ 3 , has often been used 
to integrate the ICM mass. However, the X-ray emis- 
sion is not always confirmed up to rsoo, furthermore, no 
information has been obtained about the ICM tempera- 
ture up to rsoo- These results do not ensure us to de- 
rive the ICM mass. Therefore, in addition, we introduced 
an alternative radius, risoo, within which the averaged 
mass density is 1 500-times as high as the critical den- 
sity. As shown in 4.4, this radius is as large as 0.3r2oo, 



used in Sanderson et al. (2003). We obtained the relation 

^1500 = 0.291 ( ikcv )° 65 ^50 Mpc from our sample clusters 
for kT > 1 keV, and then calculated M gas (fi5oo) as 

M gas (fi5oo) = J 4:irr 2 p dr = J 4irr 2 p c m p y 1 dr 
for each of our sample clusters. 




Fig. 11. Comparison of the gas mass of our results 
with those in Mohr, Mathiesen, and Evrard (1999) and 
Sanderson et al. (2003), for the solid square and open 
triangle, respectively. The integration radius was cho- 
sen to be the value in these two references. 



4. Correlation Studies 

LT Relation 

First, we consider the temperature and X-ray luminos- 
ity relation, which is a basic relation on clusters. In figure 
12 (left), we plot the bolometric luminosity against the gas 
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temperature. Since bright distant clusters were preferably 
proposed to be observed, luminous clusters are dominant 
in number. Nevertheless, there arc enough samples of low- 
temperature objects continuously down to kT =0.4 keV. 
We can find no clear difference in the LT relation between 
nearby (redshift z is < 0.2) and distant (z > 0.2) clusters, 
as claimed by Mushotzky and Scharf (1997), who analyzed 
the ASCA data. In the following, we consider the LT rela- 
tion by excluding objects whose accuracy of temperature 
is > 30% and redshift is > 0.5. When there are objects 
whose luminosity is by a factor of < 0.2 or > 5 higher than 
the obtained relation, we exclude them and again obtain 
the relation. 

First, in order to obtain the overall LT relation of galaxy 
clusters and comapre it with the past results, we fit the 
data at temperatures of 1.5-15 keV, by the WLS method 
(Akritas, Bcrshady 1996). As a result, we obtained the 



difference of steepness, and especially a large scatter of 
luminosity around kT ~3-4 keV, it is suggested that two 
distinct classes of objects exist; one belongs to the same 
class as richer clusters, and the other to poorer clusters; 
two classes coexist around kT =3-4 keV where higher lu- 
minosity clusters are at the cooler end of rich clusters and 
lower ones are at the hotter end of poor clusters. Such a 
break-like feature can be seen in the LT relation of Xue 
and Wu (2000). 



10° r 



10° r 



\3.17±0.15 



erg s 1 , wherc^ 10 4 
This is somewhat Ef 



relation L x = 5.80+^g x 10 42 (fcT) 
the temperature, kT, is in units of keV. 

steeper than the index of 2.72 ±0.05 reported by Wu, XueSL_ 00 . 
and Fang (1999), and near to 3.03 ± 0.06 in Xue and Wu2 
(2000), who include both rich clusters and galaxy groups. J\ 0Q 
When we chose only objects, whose luminosity is in a sim- 
ilar range of Wu, Xue, and Fang (1999), > 10 44 erg s _1 , 
the index becomes 2.61 ± 0.20, close to the value of Wu, 
Xue, and Fang (1999). In our plots, the index below 1.5 
keV becomes 4.03 ±1.07, consistent with 5.30 ±1.79 in 
Xue and Wu (2000) and 4.9 ±0.8 in Helsdon and Ponman 
(2000). Therefore, the overall relation of our results fol- 
lows the previous results, in terms of the steepness of the 
relation and one break point around kT^ 1 keV. Note that 
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X-ray faint objects around 10 



40-41.5 



erg s and fcT^0.4- 



0.8 keV arc mainly elliptical galaxies in our sample, while 
they arc galaxy groups in the previous studies. This indi- 
cates that elliptical galaxies and galaxy groups have the 
same LT relation. .2 

■43 

There are, however, several clusters whose luminosities 2 
are significantly below the typical relation around kT ^2-^ 
3 keV (NGC 4756, NGC 3258, RXJ 1833.6±652), 6-7 keV 
(A 2556B), and 10-13 keV (A 2556A); their luminosities 
are 5-10 times lower than the typical relation. The X-ray 
emission of A 2556 is doubly peaked, and the merging is 
now progressing, indicating that the temperature becomes 
higher by shock. The reason for the low luminosity of poor 
clusters NGC 4756 and NGC 3258 is not obvious; they 
might belong to the class objects on the extraporation of 
X-ray faint galaxy groups. The results of RXJ 1833.6+652 
might be contaminated by the background distant cluster 
associated 3C 383 (z=0.161; NED). On the contrary, one 
cluster, A 1885 (kT =2.3 keV), is especially highly lumi- 
nous for its temperature by a factor of 5-10. 

On the other hand, when looking at the plot in detail, 
we can see a hint of another break point of the LT rela- 
tion around 4 keV. Such a break can be seen more clearly 
in figure 12 (right), where we plot the luminosity ratio 
against the best-fit single power-law relation. When we 
obtained the steepness index separately above and below 
4.5 keV, it became 2.34±0.29 for 4 < kT < 15 keV and 
3.74±0.32 for 1.5 < kT < 5 keV. Considering this large 
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Fig. 12. Left panel is a plot of the X-ray bolometric lumi- 
nosity against the gas temperature. The three solid lines rep- 
resent the best-fit relations in the energy range of 0.3—1.5, 
1.5—5, and 4—15 keV, and the dashed line is that in 0.3—15 
keV. See the text for the equations. The right panel is the 
luminosity ratio against the best-fit single power-law relation 
(L = 5.8xl0 42 T 3 - 17 ). 



4-2. X-ray Distribution and Temperature 

The temperature dependence of the X-ray distribution 
of hot gas is an important property to consider the evo- 
lution of clusters and galaxies, since lower temperature 
systems exhibit a variety of characteristics, and it seems 
that they cannot be explained by a simple picture. We 
could derive a key issue when we understand how they 



No. ] 

connect to rich clusters concerning the X-ray properties. 
It has been reported that the lower temperature systems 
exhibit a smaller value of /? and core radius in the /?-modcl 
fitting (Mohr, Mathicscn, and Evrard 1999; Sanderson 
et al. 2003); this leads to the claim of an entropy floor 
(Ponman et al. 1999). Matsushita (1997, 2001), Beuing 
et al. (1999), Brown and Bregman (2000) indicated two 
different classes of elliptical galaxies in terms of the X- 
ray size. The ICM-to-stellar mass ratio is reported to 
be smaller for these systems than rich clusters. In or- 
der to understand these properties in a unified picture, 
we investigated the temperature dependence of the X-ray 
distribution from elliptical galaxies to rich clusters simul- 
taneously. 

In figure 13 (left), we plot the results of a single /?- 
model fitting against the gas temperature. It can be seen 
that the /? and core radius strongly correlate with the 
temperature. This is in good agreement with the result 
of Horner, Mushotzky, and Scharf (1999) and Sanderson 
et al. (2003). On the contrary, the correlation between 
the central electron density and temperature is not clear. 
Probably results of single-/? model fittings could be af- 
fected by the double structure of the X-ray surface bright- 
ness profile which is frequently seen for galaxy groups 
(Ikebe et al. 1996; Mulchaey, Zabludoff 1998); when the 
surface brightness of the double-/? structure is fitted with 
the single-/? model, the parameters /? and core radius be- 
comes smaller. In fact, as show in figure 6, some poor 
clusters apparently exhibit a double-/? model structure. 
The double-/? model is physically reasonable, since the 
depth of the gravitational potential of the central galaxy 
is comparable to that of poor clusters, and thus the hot 
gas bound by the central galaxy becomes relatively dom- 
inant at the cluster center. In figure 13 (right), we show 
plots in which we replace the results of a single-/? model 
fitting by those of outer component in the double-/? model 
fitting for the clusters listed in table . Compared with 
figure 13 (left), the correlation between /? and the core 
radius with the tempcraurc becomes somewhat small, but 
still remains. This trend is also reported by Hclsdon and 
Ponman (2000), who used the double-/? model to fit the 
surface brightness. The most remarkable difference is for 
the central electron density; the correlation with the tem- 
perature becomes clearer in figure 13 (right), in such a 
way that the correlation coefficient is 0.077 and 0.431 for 
plots of the single and double-/? model fits, respectively. In 
other words, the lower is the temperature, the lower is the 
central electron density. At the lowest temperature, some 
objects do not show the double-/? structure. The /? value 
is almost the same as that of groups and clusters, while 
the core radius is by an order of magnitude smaller. The 
central electron density is significantly higher than that 
of galaxy groups. These objects are X-ray faint elliptical 
galaxies. The difference of the X-ray surface profile be- 
tween galaxies and clusters is also indicated by Sanderson 
et al. (2003). These results indicate that we see two dis- 
tinct types of hot gas, galaxy and cluster components. 
The cluster component is dominant in rich clusters, and 
not seen in X-ray faint elliptical galaxies. The galaxy com- 
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ponent is seen in elliptical galaxies and poor clusters, and 
not seen in rich clusters. In addition, a weak correlation 
between the core radii and the temperature shown in fig- 
ure 13 (left) may be possibly due to the hot gas in the 
central elliptical galaxy. 

In summary, the X-ray distribution of hot gas depends 
on the temperature; the /? and core radius are smaller 
for lower temperature clusters, even when considering the 
double-/? model structure. This trend agrees with the pre- 
vious studies, based on the ROSAT data. On the other 
hand, in considering the central electron density, we must 
take into account the gas component, which is bound by 
the central galaxy. When looking at only the cluster com- 
ponent, the central electron density is smaller for lower 
temperature clusters. 



Figures are shown at the end of this text. 

Fig. 13. Gas-temperature dependence of the /3-model pa- 
rameters. The top, middle, and bottom panels show /3, the 
core radius, and the central electron density, respectively. The 
left panels plot the results of single-/? model fittings. The 
right panels plot the results of double-/? model fittings for the 
objects listed in table (open circle), together with those of 
single- /3 model fittings for other objects (filled square). 

4-3. The Maximum Detection Radius of X-ray Emission 

So far, the X-ray size and temperature (ST) relation has 
been reported by Mohr and Evrard (1997); the maximum 
X-ray size is proportional to the temperature. This study 
was mainly performed for higher temperature clusters. On 
the other hand, Matsushita (1997, 2001) claimed a large 
difference of X-ray size between X-ray bright and faint 
elliptical galaxies. In addition, Ponman et al. (1996), 
Helsdon and Ponman (2000), and Mulchaey (2000) in- 
sisted on the necessity of considering the undetected X- 
ray emission for X-ray faint galaxy groups. We thus have 
great interest in the connection between these classes of 
objects. The maximum radius is, of course, dependent on 
the sensitivity of the detectors, but our systematic study 
on the X-ray size of hot gas with the same instrument 
therefore has a great advantage of avoiding this effect. 
Strictly speaking, the estimation of the maximum radius 
is also dependent on the exposure time. The exposure 
time is typically in the range of 15-60 ks, and such a dif- 
ference introduces an error difference by a factor of ~2, 
and thus the difference of the maximum radius by a fac- 
tor of ~ 1.4 for the /?-model profile with /? = 0.5, when 
ignoring the background error. Therefore, considering the 
background error, the difference is smaller than 1.4. As 
can be seen from the later, such a difference is not signif- 
icant in our studies. What we have to pay attention to 
is an observed energy band, because of a wide range of 
object temperature. 

The maximum radius of X-ray size corresponds to 
-Riimit, as treated in subsection 3.5. In figure 14, we plot 
them against the temperature. Since the X-ray emission 
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of nearby objects is often larger than the GIS field of view, 
we denote such objects by circle symbols as a lower limit. 
For rich clusters, the size is as large as 3-5/i^q 1 Mpc, while 
low-temperature objects have a size of only 50-300^0 10 * 
Mpc. On the whole, the X-ray size, -Riimit , correlates with 
the temperature, as reported by Mohr and Evrard (1997). 
In addition, we can clearly see a novel feature: two breaks 
at around kT =1 and fcT^3-4 keV. Such breaks are sim-jj- 1 k 
ilar to those seen in the LT relation, and i?n m it scatters^ 
largely below 2 keV, in contrast with high-temperature ^ 
objects. Therefore, the LT relation is thought to grcatlyj 
concern the maximum X-ray size. 

The size of galaxy clusters is often defined as the radius, 
r 500 ■ although the X-ray emission has not been confirmed 
up to this radius for low-temperature objects. Then, 
we compare the X-ray maximum radius i?i; m it with ?T5oo o.Ol 
which is introduced in subsection 3.7. As can be seen in 
figure 15, the relative size of i?u m it against r^oo is quite 
different for various temperature. For high-temperature 
objects whose jtsoo is larger, i?ii m it is ~ 3 times larger 
than ri5oo- On the other hand, i?n m it is comparable to, or 
much smaller than, ?T5oo for lower-temperature objects. 
These trends are in good agreement with the comparison 
between r 500 and ROSAT R hmit (Mulchaey 2000). The 
small i?ii m it of low-temperature objects is not due to the 
low detection efficiency of the GIS for the soft X-ray emis- 
sion; the GIS count rate for the same emission measure 
is smaller by a factor of at most 2 for kT =0.5 keV, com- 
pared to kT =1 keV, and it corresponds to a reduction by 
1.4 for the detection radius for the count rate radial pro- 
file of r~ 2 , assuming (3 = 0.5. Since risoo correlates well 
with the temperature (subsection 3.7), this trend can be 
indicated from figure 14. This is strong evidence that the 
X-ray properties of galaxy clusters do not simply follow 
the scaling law, and the X-ray size of low-temperature ob- 
jects is much smaller than the theoretical prediction. We 
note that X-ray faint elliptical galaxies smoothly connect 
with galaxy clusters without a significant transition, indi- 
cating the possibility that the X-ray properties of X-ray 
faint ellipitcal galaxies and rich clusters can be explained 
within the same framework by a unified picture. 

4-4- The Gas Mass and Temperature Relation 

In treating the gas mass, we must define the integra- 
tion radius, which should depend on the system temper- 
ature. Since we have no obvious definition of the inte- 
gration radius, as described in the previous section, we 
investigated the gas mass within the two radii, i?n m it and 
r i500- For ri5oo, we did not apply a measured value of each 
object, but a calculated value from the observed relation 
of ri5oo = 0.358 (fcT/lkeV) ' 55 Mpc, so as to avoid any 
error concerned with the observed 7T500 . This integration 
radius is as large as that used by Sanderson et al. (2003), 
0.3r 200 = 0.29(fcT/lkeV) 5 Mpc. 

First, figure 16 (left) is a plot of the gas mass, 
Mgas(ri5oo), against the temperature. The M gas (ri5oo) 
correlates well with the temperature, and smoothly con- 
nects as a line from rich clusters to galaxy groups. X- 
ray faint elliptical galaxies seem to locate slightly be- 
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Fig. 14. Maximum detection radius, i?ii m it, against 
the gas temperature. The triangles represent the lower 
limit for objects whose emission is extended beyond the 
GIS field of view. 
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Fig. 15. Comparison of -Ru m it and risoo- The trian- 
gles represent the lower limit of i?n m it for objects whose 
emission is extended beyond the GIS field of view. 
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low this relation by a factor of 2-3, but it is ambigu- classes of hot gas systems appear; the low-central- mass ob- 
ous. Above kT > 1 keV, the relation is represented by ject corresponds to the class below 3-4 keV with a steep 

71 n I ~ \ /'i n _i_ n 1 ^ w 1 n!2 fi~n~* /1 \-~\ 7-\2.33±0.07 rj^i • t t — — „."u:i~ I,:™.!, — — a — ~i , , ~. i~ 



M gas {r 1500 ) = (1.0 ±0.1) x 10 12 (fcT/lkcV)" JIUU '. This 
indicates that one parameter, such as the temperature, 
determines the gas mass from rich clusters to elliptical 
galaxies. The slope of 2.33 is almost similar to the value 
1.98±0.18 for the gas mass within r 500 (Mohr et al. 1999). 

However, we should not accept the above relation 
straightforwardly, since we cannot confirm that the X- 
ray hot gas is extended up to risoo in groups and galax- 
ies, as shown in the previous section. When we replace 
A^gas (71500) by Mg as (i?i imi t) for objects whose i?ii mit is 



smaller than risoo, the relation becomes as figure 16^> 



LT relation, while the high-ccntral-mass one cprresponds 
to rich clusters with a flatter LT relation. 
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(right). The gas mass of low-temperature objects dra- 
matically decreases, and thus the relation has a break at 2 "' 
around 1 keV. Another break at around 3-4 keV is seen N - 
as the LT relation, although this also appears in figure 16" 
(left). Considering figure 16, it is claimed that the gas con- 
tent within the viriral radius is well scaled by the system 
temperature, while there is a difference in the gas distri- 
bution between lower and higher temperature objects. 

We are also interested in the scatter of the gas mass frac- 
tion. In figure 17, we plot the gas mass fraction against the 
temperature. For a higher temperature of kT > 4 keV, the 
fraction narrowly distributes in the range of 0.1-0.2h^' 5 . 
As claimed by Loewenstein and Mushotzky (1996), a sig- 
nificant scatter of the gas mass fraction within a specific 
radius indicates that the evolution of galaxy clusters is 
governed by not only gravity, but also other mechanisms. 
The fraction slightly, but significantly, decreases toward-^ 
the lower temperature, and it has a break at around kT ~ l»^o i8 
keV. This sharp break is enhanced evidence of non-gravity"-' 
effects on the evolution. I 

Last, let us look into the relation between the central 
gas density and the overall gas mass. The central gas den- 
sity is an indicator of "cooling flow rate" , and is thought 
to be related to the cluster evolution history. Here, we 
took the gas mass within 120/i^q 1 kpc of the cluster center 
instead of the central gas density, because of large un- 
certainties for the central electron densities. This is why 
we chose a radius of 120h^ kpc; this value is a median 
of the core radius in rich cluster with kT > 4 keV [figure 
13 (left middle)]. In figure 18 (left), the gas mass within 
120/i^q 1 kpc is plotted. The gas mass is not constant, but 
widely distributed by an order of magnitude. There is a 
systematic difference between lower and higher temper- 
atures; this feature should be related to that in the LT 
relation. We classify our sample by this central gas mass, 
below and above 5 x 10 n M Q , and show the ratio of the gas 
mass within risoo to the cluster-averaged one calculated 
from the relation of 1.0 x 10 12 (/cT/lkeV) 2 ' 33 " , in figure 
18 (right). Most of low-temperature objects belong to 
the low-central-density class. Interestingly, objects with 
a lower central gas mass exhibit a lower gas mass within 
the virial radius at kT ~3-5 keV, indicating that the gas 
density of a lower-central-mass object is lower at all radii 
than a highcr-central-mass object. On the other hand, no 
difference is seen above kT > 5 keV between both types of 
objects. In analogy with the case of the LT relation, two 
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Fig. 16. Gas mass against the gas temperature. Left plots 
the mass within risoo, and right plots one within ilihnlt f° r 
objects whose -Ru m it is smaller than risoo, together with the 
mass within risoo f° r other objects. The symbols are the 
same as those in figure 12. 



5. Discussion 

We analyzed about 300 objects of galaxy clusters, 
galaxy groups, and elliptical galaxies systematically with 
the ASCA data, especially the GIS data. Our results are 
free from systematic uncertainties, such as differences be- 
tween instruments and analysis procedures. In addition, 
our sample includes a wide range of system temperatures 
from 0.4 keV to ~ 12 keV. Utilizing this sample, we per- 
formed correlation studies, and found several interesting 
phenomena. Here, we discuss them from the view point 
of evolution, hierarchical structure, and so on. 
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Fig. 17. Gas mass fraction against the gas temperature 
within risoo- For objects whose -Rn m it is less than risoo, 
the fraction within -Rii m jt is plotted. 

5.1. Two Breaks in the LT Relation 

We found that the LT relation exhibits two breaks at 
around 1 keV and 3-4 keV. These two breaks are also indi- 

o 

catcd by the temperature dependence of other properties, xj 
such as the hot-gas mass, size, and gas fraction. Here, we 1 
discuss the origin and implication of the two breaks in thc^- 
temperature dependences. 

The break at around kT ~ 1 keV was already reported 
by Ponman et al. (1996), and so on, as well as the steep 
LT relation below kT = 1 keV. Ponman et al. (1996) and 
Hclsdon and Ponman (2000) suggested that this steep re- 
lation is due to galaxy feedback; lower-temperature sys- 
tems were more strongly affected by non-gravity heat- 
ing, and the gas has escaped or become thinner. Balogh, 
Babul, and Patton (1999) reproduced this steep relation 
in semi-analytical models. Alternatively, as described in 
the previous subsection, it is suggested that the X-ray 
emission of low-temperature objects is around the detec- 
tion threshold, and is possibly lost at the outermost region 
for low-surface-brightness objects, and thus the X-ray lu- 
minosity may be underestimated (Mulchaey et al. 2000). 
Accordingly, the slope of the LT relation in lower tem- 
perature objects may possibly become as small as that 
of galaxy clusters. In fact, if we estimate the X-ray lu- 
minosity within the radius risoo, rather than the detec- 
tion radius, i?n m it, the luminosity of low-temperature ob- 
jects significantly increases, while the luminosity of high- 
temperature objects slightly decreases. This behavior can 
be understood from figure 15. The slope of the LT rela- 
tion becomes 2.56 ±0.33, 3.68 ±0.37, and 3.40 ±0.59, for 
ICM temperatures of > 4 keV, 1.5-5 keV, and < 1.5 keV, 
respectively. Thus, the break at fcT~l keV disappears. If 
this is true, the existence of undetected group-scale X-ray 
emission is necessary around X-ray faint elliptical galaxies; 
we discuss this possibility in the next subsection. We in- 
fer that the hot gas escaped from low-temperature objects 
and some fractions remain, but cannot be easily detected. 
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Fig. 18. Gas mass against the gas temperature. The left 
panel plots the gas mass within 120/iJq 1 kpc. The symbols 
are the same as those in figure 12. The right panel plots 
the ratio of the gas mass of individual objects within risoo 
to the cluster-averaged one calculated from the relation of 
1.0 X 10 12 (fcT/lkeV) 2 ' 33 , by denoting objects whose gas mass 
within 120/ir. 1 kpc is less than 5 X 10 1 1 Mq as open circles. 



The second break at around kT ~3-4 keV is inferred 
from the other ICM features. The entropy and temper- 
ature relation shown by Ponman, Cannon, and Navarro 
(1999) exhibits an excess entropy below 3-4 keV. The Si- 
to-Fe abundance ratio of the ICM is constant above 4 keV, 
while it decreases toward the lower temperature below 4 
keV (Fukazawa et al. 1998). The iron-to-stellar mass ratio 
also follows the same trend (Fukazawa 1997; Finoguenov 
et al. 2000, 2001), reflecting the trend of the gas mass 
fraction that we present in this paper (figure 17). These 
three issues can be explained by the scenario that a signif- 
icant amount of ICM was pushed out by the galactic wind 
at the early galaxy formation epoch, driven by the burst 
of a type-II supernovae. Accordingly, it is reasonable that 
the break point around kT =3-4 keV is seen in the LT 
relation. The slope of the LT relation above kT > 4 keV 
is 2.34±0.29, close to 2.0 predicted by the theoretical seal- 
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ing law, while that for objects with 1.5 < kT < 4 kcV is 
3.74±0.10, significantly steep. Furthermore, a large scat- 
ter around kT ~3-4 keV implies the coexistence of ob- 
jects belonging to the higher-temperature classes in the 
flat LT relation and to the lower temperature in the steep 
LT relation. Considering these features, it is suggested 
that higher temperature objects are well-relaxed systems, 
and their non-gravity effect becomes negligible due to the 
large amount of gravitational energy. As a result, they 
follow a scaling law in the LT relation. On the other 
hand, lower-temperature objects have significantly suf- 
ferred non-gravity heating, such as galactic winds, and 
thus do not follow the scaling law. This scenario also ex- 
plains the fact that there is a correlation of the gas mass 
between the center and outer regions (figure 16). The 
non-gravity effect prevented the gas from condensing in 
the potential, and, as a result, the gas density is lower 
over the cluster-scale. 

In the bottom-up scenario, rich clusters with higher 
temperature are formed through mergers of poorer clus- 
ters. Such poorer objects are thought to belong to the 
lower-temperature classes with excess entropy. After 
merging, the excess entropy may be removed in the re- 
lease of extra-gravitational energy by radiation and so on 
to achieve virial equilibrium. Cavaliere, Menci, and Tozzi 
(1997, 1999) presented the simulation and analytic cal- 
culation of the cluster evolution, where they considered 
the preheating, subsequent hierarchical merging, shock 
heating/compression, and a new hydrostatic equilibrium. 
They found that the slope of the LT relation is 5 for galaxy 
groups, 3 for rich clusters, and saturates toward 2 for the 
highest-temperature clusters. This trend is in good agree- 
ment with our results, and indicates that smaller systems 
suffer the merger effect, while larger systems do not, be- 
cause of few merger events where comparable subclusters 
mix with eath other. 

It has been said that the LT relation of clusters of galax- 
ies is steeper than that predicted by the scaling theory, 
and thus several alternative scenarios are proposed, which 
take into account gravitational heating, radiative cooling, 
and so on (Muanwong ct al. 2002; Tornatore et al. 2003). 
Most of these scenarios predict a steep LT relation over 
any cluster temperature. However, our results show that 
rich clusters with kT > 4 keV at least satisfy the scaling 
relation. As shown in subsection 4.2, the hot-gas distri- 
bution of cluster-scale components is less concentrated for 
lower-temperature objects. Heating or cooling is a possi- 
ble effect to explain this trend. Since not all of the rich 
clusters exhibit evidence of radiative cooling, we think 
that the heating effect is preferable. In order to repro- 
duce our results, heating effects should not be so large as 
to have influence on the LT relation above kT > 4 keV. 

5.2. Gas Distribution and Two Component Models 

Concerning the X-ray surface brightness, we found that 
many objects exhibit a double-/? model structure, espe- 
cially for lower-temperature objects, while rich clusters 
and X-ray faint elliptical galaxies do not. Such a double- 
/? structure has already been found by Ikebe et el. (1996) 



and Mulchacy and Zabludoff (1998) for galaxy groups, and 
by Matsushita et al. (1998) for the X-ray bright ellipti- 
cal galaxy NGC 4636. Even for rich clusters, Ikebe et al. 
(1997) and Xu ct al. (1998) have suggested a double-/? 
structure, although the poor angular resolution of ASCA 
cannot distinguish between a double-/? model and a NFW- 
like cusp model. Following the suggestion of Ikebe et al. 
(1996) and Matsushita et al. (1998) together, it can be 
said that X-ray bright galaxy groups and X-ray bright el- 
liptical galaxies consist of the same class of objects, which 
exhibit a double-/? model structure. 

In figure 19, we plot the luminosity of the inner /3-modcl 
component against the temperature, together with the to- 
tal luminosity. It can be seen that the inner components 
and X-ray faint elliptical galaxies connect smoothly with 
each other, indicating that we might see the same type 
of X-ray hot gas. Accordingly, as described in subsection 
4.4, two types of X-ray hot gas come out; two components 
of the double-/? model are associated with the potentials 
of the galaxy and cluster. The cluster component is dom- 
inant in rich clusters, and not seen in X-ray faint elliptical 
galaxies. The galaxy component is seen in elliptical galax- 
ies and poor clusters, and not seen in rich clusters. The 
reason for this phenomenon might be due to the temper- 
ature dependence of the X-ray distribution of the cluster 
component. We have already showed that the central elec- 
tron density of the outer component in the double-/? model 
is lower for lower-temperature objects. In other words, 
the X-ray surface brightness of the outer component be- 
comes fainter for lower-temperature objects. Therefore, in 
lower-temperature objects, the cluster component cannot 
be detected and we can only observe the galaxy compo- 
nent, while in rich clusters, where hot gas of the cluster 
component is bright and the gravitational potential of the 
galaxy is relatively negligible, the hot gas of the galaxy 
component cannot be resolved. Sanderson et al. (2003) 
claimed that the properties of the X-ray halo of elliptical 
galaxies in their sample are different from those of galaxy 
clusters, and we also confirmed this trend for a larger sam- 
ple. They attributed this difference to the earlier forma- 
tion epoch than that of rich clusters. However, we suggest 
that the X-ray size of galaxy-scale hot gas is intrinsically 
compact due to the smaller scale of galaxy potential than 
that of the cluster-scale one. 

Therefore, it is speculated that, even in X-ray faint el- 
liptical galaxies, the group-scale hot gas exists, but it is 
too faint to detect. If this hypothesis is correct, we can 
claim the general view of elliptical galaxies; most ellipti- 
cal galaxies associate the group-scale gravitational poten- 
tial and thus hot gas, but the density of hot gas scatters 
widely. Figure 20 supports the above suggestion, indicat- 
ing that, for objects with kT <1 keV, the hot-gas density 
at the maximum detection radius in the GIS data is higher 
than 10~ 4 cm~ 3 , and different by only several times from 
the central hot-gas density (figure 13). In other words, 
the X-ray surface brightness of extended group-scale hot 
gas for lower-temperature objects is around the detection 
threshold, and sometimes cannot be detected even if it 
exists. In this case, we cannot distinguish whether the 
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extended group-scale hot gas exists or not by the present 

data, and Astro-E2 XIS with large effective area and low 

background level will provide an opportunity to do it. 

Our sample lacks X-ray faint galaxy groups, and HCG 

68 is the only object. The X-ray surface brightness of 

HCG 68 is centered on the elliptical galaxy NGC 5353, 

and can be fitted by a single-/? model. The maximam 

radius is at most 70/i^q 1 kpc, similar to that of X-ray faint^ 1Q 4 

elliptical galaxies. Therefore, at least the X-ray emission is u 

not due to the group scale, but due to the galaxy scale foSjooo 

HCG 68. In ROSAT PSPC observations, the extent was£. 

reported to be ~ 130/igQ 1 kpc (Pildis et al. 1995), and thus j^go r 

group-scale X-ray hot gas is indicated. In any case, the 

X-ray emission of HCG 68 seems to be dominated by the iq l 

galaxy component. Mulchaey (2000) indicated that, for 

the lower temperature galaxy groups, the X-ray emission 

becomes irregular and is dominated by individual member 

galaxies. The opposite cases are spiral-dominant galaxy 

groups, such as HCG 92 (Stephan's Quintet) (Sulentic et 

al. 1995; Awaki et al. 1997), HCG 57 (Fukazawa et al. 

2002), and HCG 16 (Besole et al. 2003). Since these 

systems do not contain a dominant elliptical galaxy, the 

X-ray hot gas of the galaxy-scale component is not bright. 

Alternatively, the X-ray emission is dominated by the faint 

diffuse component with a size of lOO-200/i,^ 1 kpc. This is 0.01 

naturally considered to be the group-scale hot gas. In 

summary, there are two types of X-ray faint objects. One 

is dominated by the elliptical galaxy in the X-ray emission, io~ 3 ^■ 

and the group-scale hot gas is hardly, or barely, confirmed""* 

The other is a spiral-dominated system, where we can 6 

observe only the group-scale hot gas with very faint diffusc^g 

X-ray emission. I 
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Fig. 19. Bolomctric luminosity against the gas tem- 
perature. The rilled squares are the same as in figure 
12, and the open circles represent the luminosity of the 
inner component in the double-/? model fitting. 



We suggest that X-ray faint galaxy groups and elliptical-^- 
galaxies may contain an amount of hot gas that is as large 
as X-ray bright ones, and there would be a non-detected 
hot gas around them. How can we observe them? This de- 
pends on the temperature of the hot gas. In our scenario, 
such hot gas would be heated up by gravity as well as 
hot gas at the inner region; therefore, the temperature is 
thought to be 0.5-1 keV. In this case, the detection of ther- 
mal X-ray emission is the most possible case, but is quite 
difficult. In the ASCA data, the systematic uncertainty of 
the CXB intensity rather than the detector intrinsic back- 
ground limits the sensitivity. Both a wide field of view and 
good imaging quality to resolve CXB into point sources 
are necessary, and such detectors are not presently being 
proposed for the future mission. The other situation is 
that the hidden hot-gas component is not virialized well, 
and the temperature is around 10 4 ~ 6 K . The detection of 
such a thermal emission is difficult, and absorption lines in 
the spectra of background objects are expected. This in- 
dication is claimed by Mulchaey et al. (1996b), and so on. 
However, we must explain why the hot gas in the outer re- 
gion is not thermalized. As described in subsection 5.1, we 
showed that the small detection radius and low luminos- 
ity of X-ray faint systems are not due to the compactness 
of the hot-gas extent, but due to the detection sensitiv- 
ity. Several numerical simulations show that the ICM is 
heated up by shock waves within a radius of rsoo (Evrard 
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Fig. 20. Electron density at the maximum detection 
radius, Rumit- The open circles represent the case 
where the emission is beyond the GIS field of view. 
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et al. 1996; Takizawa, Mineshige 1998). Therefore, it is 
reasonable to think that the hot-gas temperature is not so 
different from the inner region as objects whose detection 
radius is large. 

5.3. Scatter of X-ray Luminosity of Galaxy Groups and 
Elliptical Galaxies 

Ao far we suggested, two hot-gas components exist in 
lower temperature systems: a galaxy component and a 
group component. Considering these components, we dis- 
cuss the scatter of the X-ray properties for elliptical galax- 
ies and galaxy groups. As mentioned in subsection 5.1, we 
imply that these two systems are basically the same class 
of objects for hot gas and dark matter, and only the stellar 
distribution is different; stars in elliptical galaxies concen- 
trate in one galaxy, while stars in galaxy groups separately 
locate in several galaxies. This difference may be due to 
whether the galaxies in the systems have merged or not. 
The scatter of the X-ray luminosity is mainly caused by 
the scatter of the hot-gas density for the group compo- 
nents and the brightness of the galaxy components. The 
origin of the former is related to the degree of heating, 
dark-matter concentration, system age, and so on, and 
may be explained in formation theories. On the other 
hand, the scatter of the X-ray brightness of the galaxy 
components is not simply understood. Matsushita (1997, 
2001) showed that the hot-gas mass within the 4-timcs 
effective radius is systematically different between X-ray 
bright ones and X-ray faint ones. Therefore, the differ- 
ence at such a small scale is thought to be attributed not 
to the group component, but to the cD elliptical galaxies. 
The compression of the galaxy hot-gas component by the 
high pressure of the group hot-gas component in X-ray 
bright ones, the escape of hot gas from X-ray faint ones, 
and so on, are considered, although we cannot investigate 
this issue with the ASCA data. High-resolution imaging 
with Chandra data will give an answer. 
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Fig. 13. Gas-temperature dependence of the /3-model parameters. The top, middle, and bottom panels show (3, the core radius, 
and the central electron density, respectively. The left panels plot the results of single-/3 model fittings. The right panels plot the 
results of double-/? model fittings for the objects listed in table (open circle), together with those of single-/3 model fittings for other 
objects (filled square). 
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